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We study two numerical techniques based on the homotopy perturba- 
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involving the Caputo—Prabhakar fractional derivatives. In this manuscript, 
the convergence for numerical solutions obtained using HPTM and the con- 
vergence and stability for numerical solutions obtained using FABM are inves- 
tigated. We compare the solutions obtained by the HPTM and the FABM for 
some nonlinear time-fractional differential equations. Moreover, some numer- 
ical examples are demonstrated in order to show the validity and reliability 
of the suggested methods. 
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1 Introduction 


In recent years, fractional calculus, which is a generalization of the integer- 
order calculus or integration of integer order to arbitrary order (real or com- 
plex order), plays an important, fundamental and significant role in mathe- 
matical modeling and solving the models obtained in different fields of math- 
ematics, sciences, physics, finance and artificial intelligence, and engineering; 
see [11, 39, 3, 10, 31, 38, 2]. Nowadays, fractional differential equations have 
been applied extensively as an important and crucial tool to describe many 
different types of complex mechanical, physical behaviors, sciences, engineer- 
ing, anomalous diffusion, vibration, control, viscoelasticity, electrochemistry, 
and others, since most of the models in different research areas and engineer- 
ing applications are nonlinear and obtaining their solutions has an important 
role. 


In this manuscript, we consider a numerical method based on the ho- 
motopy perturbation transform method (HPTM), which is a combination of 
homotopy analysis method and Laplace transform scheme [29, 40] for solving 
the following nonlinear time-fractional partial differential equation (FPDE) 
with m-—1<f<m: 


“Dy wor U(t) + Au(t) + Nu(t) = A(t), 
u)(0) = uk, k=0,1,2,...,m—1,meEN, (1) 


where Au(t) and Nu(t) are linear and nonlinear eperavor applied to u(t), 
respectively, and h(t) is a continuous function. Also © 4 uw,0+ denotes the 
Caputo—Prabhakar fractional derivative defined by Garra et al. (see [11]) 


o Duis u(t) = i: (t— 7) # Eilat — 7)? )t(r)dr 


= d 
=E,T u,w,0+ Gp u(t), O<p<l, (2) 


where E is the Prabhakar. It is obtained by modifying the Riemann— 
Liouville integral operator [31] by extending its kernel with a nonsingular 
three-parameter Mittag—Leffler function (ML), which is known as the Prab- 
hakar function and its integral is (see [11]) 


(EY w(t) = | (t— 7)! BY, (w(t — 7) )u(r)dr, 


where £7, is the Prabhakar function and its function is (see [11]) 


Bt io= >, et Rp), Ku) >0,y>0, tec, (3) 
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and (7),, is the Pochhammer symbol, which is formulated as [15] 


(Yo=h ), =774+1).--(¥+n-1), neEN. 


Most of our interest in studying ML functions is related to their importance 
in fractional calculus which arise in modeling a complex susceptibility in the 
response of disordered materials and heterogeneous systems [26], the response 
in anomalous dielectrics of Havriliak-Negami type [12, 37], in fractional vis- 
coelasticity [14], in the discussion of stochastic processes [9], in probability 
theory [16], in the description of dynamical models of spherical stellar sys- 
tems [1], and fractional or integral differential equations [25, 6, 5, 7, 8, 4, 18]. 
In this manuscript, we consider the following nonlinear time-fractional order 
differential equations with initial conditions: 


. DP wwo,0+ Ut) = f(t), 
u*“)(0) = uk, k=0,1,2,...,m—1,meEN, (4) 


and in order to solve the nonlinear FPDE (4) with m—-—1 < wp < m, 
we apply a numerical technique based on the fractional Adams—Bashforth 
method (FABM) to obtain the numerical solution of (4). The present work 
aim is to compare the numerical solutions that are obtained by applying 
the FABM and HPTM to solve nonlinear time-fractional order differential 
equations. Several numerical-analytical methods have been applied to solve 
linear and nonlinear differential equations of fractional derivatives in order 
to obtain exact solutions and numerical solutions; for instance, the Her- 
mite collocation method [30], optimal homotopy asymptotic method [17], 
g-homotopy analysis transform technique [36], homotopy analysis Laplace 
transform method [23, 24], the homotopy analysis Sumudu transform method 
for solving the linear and nonlinear Fokker-Planck equations [33], and other 
methods [19, 20, 28, 27, 34, 21] were suggested. 

This paper is organized as follows. Some necessary definitions and mathe- 
matical preliminaries of the fractional calculus are introduced in Section 2. In 
Section 3, we introduce a numerical method based on the HPTM. In Section 
4, the FABM for time-fractional linear and nonlinear differential equations 
is introduced to obtain the numerical solution. In Section 5, the proposed 
methods are applied to some nonlinear time-fractional differential equations 
to verify the validity and applicability of the proposed methods. Finally, a 
brief conclusion is added at the end of the manuscript. 


2 Preliminaries 


In this section, we study some important and basic properties of fractional 
calculus theory such as Laplace transform, definitions, and lemmas, which 
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are applied in the next sections. 
Definition 1. [see [22, 31]] Let 0 < a < 1, let f € L'{a,0], and let 


0<t<b<o. Then the left-sided and the right-sided Riemann—Liouvillel 
fractional integrals and derivatives of order a are defined as follows, 


I, f(t) = _ / f(r)(t—1)°Nar, 


b 
Ie f(t) = _ | flr)(r —t)°*ar, 


Daf) = TG f LOe— rar 
b 
DMO) = Ty f fOr 


where ['(-) is the Gamma function. Also, the left-sided and the right-sided 
Caputo fractional derivatives of order a are given by 


7 oe eee : sae 
Dal = Ne Gl = pqaw f tee senar 


b 
ODE f(t) = —Ho* SF) = aren kd ie fr)ar. 


Lemma 1. [4] Let p, u,w,y € C, let R(u) > 0, and let R(p) > 0. Then the 
Laplace transform of (2) for m—1< ps < m is given, 


m—1 

c(CD u(t); s) = s*(l1—ws ?)7U(s) — sh—F-1 (1 — we?) 1u* (0+), 
Pb, ,0 = 

(5) 


where U(s) = [,° e~*‘u(t)dt. Also, the Laplace transformation of the Prab- 
hakar function (3) is denoted by (see [11, 15]) 


L(t Bt, (wt”); s) =s *(l—ws-?)7”, 
Lemma 2. [22] Let p,u,w,y € C and let R(w) > 0, R(p) > 0. Then, 
t 
/ (f= y)¥ "ER (w(t — y)?)y’ "dy = Pvt?" BF ny (wt?). (6) 
0 


Lemma 3. [22] Let p,p,v,w,0,7 € C and let R(u) > 0, R(p) > 0, R(v) > 0. 
Then the following relation is hold for any summable function u € L(a, b), 
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a — pyto 
51,00 Ep,y,0,0+U ~~ Be ptt u (7) 


EE’ 
p 
Also, by substituting 0 = —7¥ in (7), the following relation is obtained, 


E’ E 


_ ppv 
Pp, H,w OF pw orl = Lor me (8) 


3 HPTM to solve the nonlinear fractional order 
differential equations 


In this section, we obtain a combined numerical method based on the ho- 
motopy analysis method and the Laplace transform method for solving the 
nonlinear FPDE in Caputo—Prabhakar fractional derivatives sense. Using 
formula (5) for m = 1 and applying the Laplace transform on both sides of 
equation (1), we get 
(ODF , «oti + Au(t) + Nu(t) — A(t); s) =0, 
s#(1 —ws-?)7U(s) — s#—1(1 — ws-?)%u(0T) + £(Au(t) + Nu(t); s) = H(s), 
1 
U(s)=s *(1—ws °)-7H(s) + —u(0) — (s “(1—ws ”)77) 
S 
x [< (Au(?) + Nu(t); s)|. (9) 
By applying the Laplace inverse transform on (9), we obtain 
1 
u(t) =£72 Faw —ws"?)7H(s) + =u(0); t ta (#1 —ws-?)-7) 
x [< (Au(?) +Nu(t); s)| ; q ; (10) 


To obtain the solution of (10), we use the HPTM. Suppose that the solution 
of (1) is u(t), which can be expressed as the following infinite series: 


u(t) =} punt), (11) 
n=0 


where u,,(t) for n = 0,1,2,... are known functions. Also, the nonlinear term 
of (1) is Nu(t), which can be presented as the following infinite series: 


Nu(t) = 5° p"H,,(u(t)), (12) 
n=0 
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where H,,(u(t)) = {4 oo [V( yk viui)| a} are polynomials of the Ado- 


mian polynomials type; see [13]. Now, substituting (11) and (12) into (10), 
we conclude that 


5 pun(t =H [(s*(1 — ws ?)-7) 
x [c(Aprunit) + No p"un(t)ss) A], (13) 
n=0 ar 


where H(t) = £7! [s-#(1 —ws’)~7H(s) + +u(0); i . By equaling the coef- 


ficients on powers of p on the both sides of relation (13), we obtain the series 
solution as follows: 


ps uo(t) = H(t), 
ph su (t) = —£-1[(s-# (1 ~ ws7P)-7) x [C(Aao(t) + Ho(u(t)):s) [se], 
ps un(t) = £71 [(s-M(1— ws ?)-7) x [£(Aur(t) + (ul): 8) |i], 


pt): ungi(t) =—L ‘(s #(1 —ws7?)-7) x [£(Aun(t) + Hn (u(t)s8) |e]. 


So, the solution of (1) by using the HPTM is obtained as follows: 
u(t) sae > ug(t) = ape mel 


Theorem 1. Let S,, = S079 u(t) and let u(t) = O72 9 ux(t) such that 
u(t), un(t) € C((0, \)) and for any n = 1,2,3,..., || Un+i(t) I< A |] un(£) II, 
0<2A<1. Then the solution of (1) using the HPTM is converges. 


Proof. To prove convergence, we need to show that the series S;, is a Cauchy 
sequence. We have 


n+1 n 
Il Snt1 — Sn |] =I 55 u(t) — $0 ue(Z) | 
k=0 k=0 


|| unt IIS A |] un) IS A? |] una) IIS 
<A"*" || uo(t) [I - 


Then for any m,n € N,n > m, we obtain 


I| Sn — Sm || 
=|| (Si, Sn—1) t (Sn-1 Sn—2) t (Sp-2 Sin—3) a i (Sin41 _ Sim) | 


Comparison of homotopy perturbation transform method and ... 69 


<|| Sn, = Sn-1 | T | Sn-1 — Sn—2 | = | Sn—2 aa Sn—3 | + | Sm-+1 =, Sm | 
<A" || wo(t) |] +A"? |] wo(E) |] 4-2. +A" |] wo (Z) || 
y—)\"-™ i \nrri 
nD WLC oe _ \yn-m : 
< "|| wold) IIs Il volt) |, 1-2" <1 
As m,n — oo then || S;, — Si, || 0, since uo(t) is bounded and 0 < A < 1. 


So the series S;, = a uz(t) is a Cauchy sequence in the Banach space and 
therefore converges. 


4 The Adams—Bashforth scheme to numerically solve 
nonlinear fractional order differential equations 


In this section, we consider the fractional differential equation (4) with initial 
conditions and apply the Adams—Bashforth scheme in order to numerically 
solve this equation. Using relation (8) in Lemma 3 and applying the inte- 
gral transform operator Be onaiine to both sides of the fractional differential 
equation (4), the solutions of this equation are obtained as follows: 


ES w,0+ (CDF 20+ ult) = ae f(t), 


By ioe an qu) = ED word (t)s 
u(t) — wo) = f (er BR, ole rf 
u(t) = u(0) + i (t— T)Y TED (w(t —7)?)f(r)dr. (14) 


Now, we first put t = ty4i, n =0,1,2,3,... and thent =t,, n =0,1,2,3,... 
in relation (14). Now we obtain 


u(tn 41) = u(0) + | ting — TEED (w(tnga — TF (nde, (18) 


tn 
tu(tn) = u(0) + | (ta — 7) LEY (lta — 7)?)F (rar. (16) 
0 
By subtracting two relations (15) and (16), we obtain 
bei 
nee ae | (tn4a — 7)" EY, (w(tnga — 1?) F dr 
0 


Z | "(tn — T)HLEY , (w(tn — 7)°) f(r). (17) 
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To solve (17), we approximate the function f(r) by the Lagrange interpola- 
tion as follows: 


T—tn-1 T—tn 


P(r) =; = 


P(r)s 


f (tn) + 


T—tn-1 T—tn 


Ht F(t) — 


Fltn—1), 


tn—-1 = tn 


f(tn-1); (18) 
tj = jh, 7 =0,1,2,3,..., h= tn — ty-1. 


Thus, by substituting (18) in (17), we obtain 


u(tn41) =ultn) His 
+ [Ae FO tan PN lags = YT ta) 
= eg PO gs 1 (elton — TP YCr = tn] 
= [AE Pt = 2 ella — 7) Yr = ta 


= ee PO tn = 1 1B alte = TPCT — tar] 


til) + [fe a Mina) 4, | as ea 


where 


tn41 
Ay = / (trai — eh A obaa — T))(7 — ty_1)dr, 21) 
0 


tn4i 
ae i (toga — TYE? (w(tn4a — 7)°)(7 — ta) dr, 22) 
0 
tn 
Ape (ty — 7)" BY, (w(ta — 7)°)(7 — ta—a)dr, 23) 
tn 
Ana | (ta — TY 1BY,, (w(ta — 7)°) (7 — ta) dr. 24) 


To calculate relations (21)—(24), first we calculate the integral (21) as follows: 
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tn41 
Ana =| (tn4i — T)Y TE? (w(tn4a —T)P)(7 — ty_-1)dr 


tn41 
-| (tn41— 1m ame OLE (Clare —7))rdr 
tn4i1 
—tn—1 | (tn41— TY BY (w(tn4a —T)P)dr. (25) 
‘ 


By applying (6) in relation (25), we obtain 


ted 
Pie | (toga — T)#TEY,, (w(tnga — 7)°)(7 — tna) 
0 
D2) ing) Ee ase aya)?) tack X PO Ga) PE) sl Gaga)”) 


=(n+ DP ne BD saan + 1)h)?) 


= (n= Ah x |(n + 1)HHED 41 (w((n + hy?) 
a(n + 1h? x [(n + YEP o(w((n+ Da)”) 
— (2 YEP p(n + Yh)? |. (26) 


Similarly, for the residual terms, we have 


tadi 
Ana =f (tnaa — 7B} leeltnss — TY)(r— ta) 
=(n + 1)tait! 


x (n+ DEP po(w((n + h)*) = EP yr (w((n + Yh)")], (27) 


tn 
An3 = / (tn — PTET Lait —7)?)(7 — ty_1)dr 
0 


= ntpett x [nz 


7 pra(w(nh)?) ~ (nL) EP vya(w(nh)®)], (28) 


p+ 
tn, 
die | (ta — #1 BY, (w(tn — 7)°)(7 — ta) 
0 
= nitlpetl y [EP +2((mh)”) = FP 41 (w(nh)?)| (29) 


Substituting relations (26)—(29) in (19), we obtain the approximation scheme 
to solve numerically (4). 


Theorem 2.[Convergence] Let u(t) be a solution of (4) and let f(t) be con- 
tinuous and bounded that is introduced in (4). Then the numerical solution 
of u(t) is given by 
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tn41 
u(tn41) =U(tn) + } (tn41 — TPT ED y(w(tn+1 =e )") 
0 


T—tn-1 oS tn 
iF = tp-1 P(tn) tn—1 = Ff (tn1)| a 
tn 

= [etn — TEL leoltn = 79?) 

0 

T= tn—1 T= tn 

Ee _ tn—-1 I(tn) | tH = Atma) ial 

+ R(t, 4,7), 


where R(t, u,7) is the error function and || R(¢, w,7) |lao< M. 


Proof. Using (17) and (18), we have 
tn41 
ti(tn41) =U(tn) + | (taga — tT)! ET, (w(tn4a — 7°) F (rar 
0 


2 i "(tn — 7) LED (w(t — 7)°)f (nar 


=u(tn) 
tn41 
+ | (tn4a — TPT EY (w(tna —r)?) 
0 
T—tn-1 T— ty fTtD(r) < 
E an tn-1 H(tn) . tn—1 _ gf tn—) : (n + 1)! I i) ar 


7 i "(tn — TH LED, (w(tn — 7”) 
0 


T—tn-1 T—tn ()(7) 
[EH ttn) + tin) + FD TT or - aa] ar 


tn = tn—1 tn—1 = tn ( )! i—0 
=u(tn) oF Lit, L, n) ae R(t, HL, n), 
where 
L(t, 1,7) 
tn41 m 
=f eat — 2S eoltnts — 7) 
T—tyh-1 T—ty 
——— f(t, ——— f(t,_1)|d 
Lan tn) + pS tea) ar 
tn 
- ‘| a= ee olay, —T)?) 
T—tyh-1 T—tny 
tn) 4 tn—1)| dT, 
Laan te) + pS tea) ar 
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R(t, 1, n) 


tet (n+1) ie 
= [nt = ER loltns — 7) (EP 


(n+ 1)! a 


(r ti))dr 


a w—-l py p FG) aT : 
= (tn — TTB (woltn = 7))( [I -#))ar. 


(oars 


To establish convergence, we need to show that || R(t, u,7) |lo< M, MEN. 
For this, we have 


| R(t, 4,7) lhoe 


tn41 (n+1) r) we 
=f naa — BY etna — (FP 


(n+1)! 2 


— [arm otter) (LL TY oe — tar 


0 (is 
tn41 (m4) (7) 
<M ts PBR elton — 9) Taye LD ~ 10) 2 Ne 
* 4=0 


+ || rf (tn — "1 BY, (w(t ny (= r) in bar i 
[reo] (t) | ~ 
“cena (vad)! | Ie i) loo 
tn+1 
I (tng — T)Y EF (w(tnga — 7)?) dr 
i | i en ie (tn —T)* 1 ED ,(w(tn — 7)?)dr 
te (0, a (n) oo n yh a ; 
Using Lemma 2, we obtain 

[70] m+n()| 


R(t oo oo thai 
| ( |, 7) | < te[Otovs] ( (n+ 1)! | He i) ll C41 Rermee thei) 


ie dl 
| Th (t— ti) lloo thE Dwi (wth) 


ma 
te(0, a 
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Fora Frc 
< su { max ——WW—, max —\ 
Bea t€[0,tr41] (N +1)! 'tefotryi] nl 
nlartt (n—1)!h" 
x 4 tip, uti (thir) - a ET a (wth)) 
= M. 


Theorem 3.[Stability condition] If 


> 0 


M, nth” i + 1) th FP ta (Wty) , thE? +1 (wth) 
4 h h2 


as n + oo and h > 0, then, the FABM is stable. 
Proof. Using (17), we have 


I| w(tnt1) — U(tn) lloo= II | (trp — TP ED y(w(tmsa — 7)?) f(r) dr 


: i "(ta — TH LED, (w(tn — 7)? (7) AF loo 


< || | , (ee = TPT EY (w(tns1 —1)?) f(r)dr loc 


+f "by — ERY, ulin — 7)? (AT lle 
Hence 
u(tne1) = (tn) [lo 
hf Guts = lw oe cele ec 


—t 
+f G9 BY Gt, —*7) (ST a ame: 
<|| P(t, 4,7) lloo + |] Rn (6) lloo: 
where || P(¢, 4,7) ||. and || R4(t) ||. are calculated as follows: 
[P(t, H,7) lle 


=I ) Cb — TTY (w(tnga — ry) II ENE AT loo 


k=0 k=0 
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te) lle tei? are ee 
| f( k) | re Dei n+ UIT | te 
k=0 
Su! (tx) lloo th n+ ED gi (wth 41) rth” 
= 4h 
and 
|Ri(t) |hoo 
SH rated flte) 
— tr k 
=I [G TYE a ett — 7Y)| 0 TD ay | ls 
a ae re 
= || FG) Ie # dal 2 eae 
<y! (tk ) ll n spei( "TT | tr| 
k=0 k=0 
aN = Il F(t) loo HHEF 41 (wth)(m — 1)!ar— 
- 4h ; 
k=0 
Then 
Iu oe 4 u(t n) lloo 
> | f(t) lloo th meee, ee (wth 4, )nth” 
4h 
k=0 
n-1 n—-1 
a5. Il Fti) loo thED 41 (wth)(n — 1th 
4h 
k=0 
M, nth” i + thar Eb wp (trg) 2 en) 
4 h h2 : 


where My = max;ejo,t 


f(t) . Therefore, 


n4+1] 


M, nth” iE mt 1) ED 41 (wt) ; nes) 


m rs he = 


asn— oo, h > 0 in which h = =O is considered. 


5 Illustrative examples 


In this section, we compare the obtained solution with the numerical sim- 
ulations of the applied methods on five test problems in order to show the 
accuracy and efficiency of the proposed techniques. 
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For the special case whenpu=0.95 
Exact 
2.57 ——HPTM 
=== FABM 


u(t) 


u(t) 


M.H. Derakhshan and A. Aminataei 


For the special case whenu=0.99 


Exact 
- - -FABM 
—— HPTM 


Figure 1: Comparison between numerical solution and the exact solution of Example 1 


when p = 0.95, 0.99, p=w=1. 


x10% 
3 
= SU aU ral for 10.95 | 
1 
. 2 7 - 
re) ’ re) 
ae _ 
= = 
Ww / Ww 
1 
oc 
oe = id 
0 1 
x 107 
6 
= = = UEyact4 Here! for #=0.99 } 
- 4 
fe} 
a 
= 
Ww 
2 ri 
’ 
’ 
oe “se 2 
0 1 3 


x10% 
wane toa We yact™ Fag! for #=0.95 F 
2 i 
0 Pal i al 
0 1 2 3 
t 
x10% 
6 3 
vero [Uevact Yeap! for #=0.99 F 
= 4 $ 
° : 
= : 
_ : 
Bi 2 
2 7 
i erases” 
0 1 2 3 
t 


Figure 2: The value of absolute error of Example 1 when p = 0.95, 0.99, p=w=1. 


Example 1. We consider the following nonlinear fractional order differential 


equation: 
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Table 1: The exact solution, the numerical solutions, and the absolute errors 
of Example 1 for different values of t when pp = 0.95, p=w=1. 


t UBeact(t) | UyprM(t) | UraBa(t) | |UBeact(t) — unPrM(t)| | Ueeace(t) — UraBM (t)| 

0 0 0 0 0 0 
0.01 .0001 0.0011 .0008 0.001009960804482 0.000677271813138 
0.02 0.0004 0.0024 0.0017 0.002022248070965 0.001297961649675 
0.03 0.0009 0.0039 0.0027 0.003027289212479 0.001857141698735 
0.04 .0015 0.0056 .0039 0.004020110685205 0.002353085479644 
0.05 0.0024 0.0074 0.0052 0.004997320737974 0.002785155766582 
0.06 .0034 0.0093 .0065 0.005956376904834 0.003153291833384 
0.07 .0046 0.0115 .0080 0.006895275398590 0.003457792029013 
0.08 0.0059 0.0137 0.0096 0.007812388169833 0.003699199718883 
0.09 0074 0.0161 .0113 0.008706366104405 0.003878235523084 


Table 2: The exact solution, the numerical solutions, and the absolute errors 
of Example 1 for different values of t when pu = 0.99, p=w = 1. 


t UBeact(t) | UyprM(t) | UraBa(t) | |UBeact(t) — er |UBeact(t) — UraBM (t)| 
0 0 0 0 0 0 
0.01 .0001 0.0003 .0002 0.000201585080800 0.000111408806560 
0.02 .0004 0.0008 .0006 0.000403349846034 0.000164732892224 
0.03 0.0009 0.0015 0.0010 0.000603578912785 0.000160544488950 
0.04 .0015 0.0023 .0016 0.000801354185865 0.000099955930106 
0.05 0.0024 0.0034 0.0024 0.000996039906343 0.000015740815560 
0.06 0.0034 0.0046 0.0032 0.001187154171698 0.000185164079811 
0.07 .0046 0.0059 .0042 0.001374313737802 0.000406881133539 
0.08 0.0059 0.0075 0.0052 0.001557204785139 0.000679428650402 
0.09 0074 0.0091 .0064 0.001735565417799 0.001001324957541 
ChyyY _ 2 
De a,o+ Ut) = f(t) —uU (t), 
subject to the initial condition 
u(0) = 0 
and 
pth 


f(t) 


u(t) =[2¢(1 — + eae 


=(#-3P° —# 


2 3t?(2¢(1 


Under these conditions, the analytical solution is given by (see [41]) 
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We apply the proposed schemes in this article and solve this problem. Figure 
1 shows the analytical and numerical results of Example 1 obtained by ap- 
plying the HPTM and the FABM together with the exact solution. Also the 
absolute error for 44 = 0.95, 0.99 with several values of t is shown in Tables 
1 and 2, and for w = 0.95, 0.99 with several values of t, it is shown in Figure 
2. The numerical results demonstrate that the FABM method is an efficient 
algorithm. 


For the special case when p=0.95 


6 T 


u(t) 


For the special case when p=0.99 


Exact 
——FABM 
——HPTM 


Figure 3: Comparison between numerical solution and the exact solution of Example 2 
when p = 0.95, 0.99, p=w=1. 


Example 2. We consider the following nonlinear fractional order differential 
equation: 


“Dwr Ut) = F(t) — u(t), 
subject to the initial condition 
u(0) =0 
and 


7 2, yaa, 26 * 67 Ge 
fG) = (te -37 +) 4 T@om Fo aT Ta>p’ 


Under these conditions, the analytical solution is given by (see [41]) 


u(t) = 2¢ -3 +2. 
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[Ue xact™UHpTm! for 4=0.95 


Error 


x10% 


Ue act UHpT! for w=0.99 


Figure 4: The value of absolute error of Example 2 when pp = 0.95, 0.99, p= w = 1. 
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For the special case when »=0.95 
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Exact 
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Figure 5: Comparison between numerical solution and the exact solution of Example 3 


when p = 0.95, 0.99, p=w=1. 


We apply the proposed schemes in this article and solve this problem. Figure 
3, shows the analytical and numerical results obtained by applying the HPTM 
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For the special case when u=0.95 


— Exact 
——HPTM 
- - -FABM 


For the special case when p=0.99 


— Exact 
——HPTM 


u(t) 


- - -FABM 


Figure 6: Comparison between numerical solution and the exact solution of Example 4 
when p = 0.95, 0.99, p=w=1. 


For the special case when p=0.95 


Figure 7: Comparison between numerical solution and the exact solution of Example 5 
when p = 0.95, 0.99, p=w=1. 
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For the special case when p=0.95 


log, 0 Error 


log, 0 Error 


Figure 8: The value of absolute error of Example 5 when ps = 0.95, 0.99, p= w = 1. 


and the FABM together with the exact solution of Example 2. Also the 
absolute error for 4 = 0.95, 0.99 with several values of t is shown in Figure 
4. The numerical results demonstrate that the FABM method is an efficient 
algorithm. 

Example 3. We consider the following equation: 


2t?-# 
Diao or elt) + ult) = 5 


— 2 = < 1. 
P;H,wW 3—p) +t ? u(0) 0, 0< HS 


For this example, the exact solution u(t) = ¢t? was presented in [32]. This 
example is solved by using numerical methods introduced in this article for 
different values of , and their diagrams are shown in Figure 5. 

Example 4. We consider the following equation: 


1 1 t 
VF pw or W(t) 5 u(t) = 50 U5): u(0) =1, 0<p<1,0<t<10. 


The exact analytical solution of this example is u(t) = e', which was presented 
in [18]. This example is solved by using numerical methods introduced in this 
article for different values of 4, and their diagrams are shown in Figure 6. 
Example 5. We consider the following nonlinear fractional initial value 
equation: 


CDT yoru(t) + u(t) =f(), u(0)=1, 0<p<i,0<t<1, 
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4032008-# = ==T(5+4)t4-2 or(u+1 3t2 3 
ft) = \ 2)t u 4 i] ; 
T(9 — p) T(5— 4) 4 


The exact solution of this example is u(t) = ¢8 — 3¢4+2 + %", which was 
given in [35]. Diagrams of the numerical solutions for different values of wu 
are shown in Figure 7. We calculated the logarithm of the absolute error for 
different values of , and their diagrams are shown in Figure 8. 


6 Conclusions 


In this paper, the FABM and the HPTM were presented to study the initial 
value problem of the linear and nonlinear fractional differential equations 
involving the Caputo—Prabhakar fractional derivatives for numerically solving 
this type of equations. In addition, the two proposed methods were compared. 
Numerical examples and graphical representations were presented for testing 
the validity and accuracy of suggested methods. These examples showed that 
the FABM is a strong tool for obtaining the numerical solutions of linear 
and nonlinear fractional differential equations. Furthermore, the stability 
property of the FABM and the convergence analysis for both methods were 
considered. The numerical results were stable and converged well for different 
values of t for all examples. 
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